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The mechanisms of thermally activated magnetization 
switching in small ferromagnetic particles driven by an ex- 
ternal magnetic field are investigated. For low uniaxial 
anisotropy the spins rotate coherently while for sufficiently 
large uniaxial anisotropy they behave Ising-like, i. e. the 
switching then is due to nucleation. The crossover from co- 
herent rotation to nucleation is studied for the classical three- 
dimensional Heisenberg model with uniaxial anisotropy by 
Monte Carlo simulations. From the temperature dependence 
of the metastable lifetime the energy barrier of a switching 
process can be determined. For the case of infinite anisotropy 
we compare numerical results from simulations of the Ising 
model with theoretical results for energy barriers for both, 
single-droplet and multi-droplet nucleation. The simulated 
barriers are in agreement with the theoretical predictions. 



With decreasing size ferromagnetic particles become 
single-domain which improves their quality for magnetic 
recording. On the other hand, when the particles are 
too small they become superparamagnetic and then due 
to thermal fluctuations no magnetic information can be 
stored (see e. g. for a review). Hence, much effort 
has been devoted to an understanding of the behavior 
of small magnetic particles experimentally , analyt- 
ically J5|, and in computer simulations. In the latter 
case, mainly magnetization switching by nucleation |]-|^] 
has been studied using Ising models of finite size, but 
also other reversal mechanisms in models with continuous 
degrees of freedom like coherent rotation, single-droplet 
nucleation, and multi-droplet nucleation have been dis- 
cussed §,§. 

In this paper we focus on different thermally activated 
reversal mechanisms occuring in ferromagnetic particles. 
We will consider a finite, spherical, three-dimensional 
system of magnetic moments. These magnetic moments 
may represent atomic spins or also block spins following 
from a coarse graining of the physical lattice Our 
system is defined by a classical Heisenberg Hamiltonian, 



H 



axial anisotropy which favors the z-axis as easy axis 
(anisotropy constant d > 0). The last sum is the cou- 
pling of the spins to an applied magnetic field, where 
B is the strength of the field times the absolute value 
of the magnetic moment of the spin. We neglect the 
dipolar interaction. Therefore, the validity of our results 
is restricted to particles which are small enough to be 
single-domain in the remanent state pd[ | . 

In the following we will investigate the thermally ac- 
tivated reversal of the magnetization of a particle which 
is destabilized by a magnetic field pointing into the di- 
rection antiparallcl to the initial magnetization which is 
parallel to the easy axis of the system initially. In this 
case after some time the particle will reverse its magne- 
tization. 

Due to the many degrees of freedom of a spin system 
numerical methods have to be used for a detailed micro- 
scopic description. Since we are especially interested in 
the thermal properties we use Monte Carlo(MC) meth- 
ods for the simulation. Although a direct mapping of the 
time scale of a MC simulation on experimental time scales 
is not possible this method provides information on the 
dynamical behavior since it solves the master equation 
for the irreversible behavior of the system. 




where the Si are three dimensional vectors of unit length. 
The first sum which represents the exchange of the 
spins is over nearest neighbors with the exchange cou- 
pling constant J. The second sum represents an uni- 



FIG. 1. Snap-shot of a simulated spin system during coher- 
ent rotation. In the picture we omit the interior part of the 
system and show only the outer spins of the sphere. R — 6 
spins, B = 0.7 J, d = 0.35J/spin, T = 0.09J. 

We consider spins on a simple cubic lattice of size 
L x L x L and simulate spherical particles with radius 
R = L/2 and open boundary conditions on this lattice 
using the MC algorithm described in B. We start our 
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simulations with an initial spin configuration where all 
spins are pointing up ( Si = (0, 0, 1)). The magnetic field 
B = (0, 0, —B) destabilizes the system and after some 
time the magnetization will reverse. The metastable life- 
time r is defined by the condition M z (r) = where M z is 
the z-component of the magnetization M = (1/N) J2i S». 

For sufficient low anisotropy the spins can be expected 
to rotate coherently. Such a reversal process is shown in 
Fig.|l| where a spin configuration of a simulated system 
of size R = 6 spins during the reversal process is repre- 
sented. The spins are nearly parallel, except of thermal 
fluctuations. Following the theory of Neel Jl2j and Brown 
plU the energy barrier which has to be overcome during 
the reversal is only due to the anisotropy of the system, 



AE r 
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(2) 
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where M is the saturation magnetization of the system 




FIG. 2. Snap-shot of a simulated spin system as in Fig. 
[j] but for single-droplet nucleation. The z— component of the 
spin is color-coded (lighter grey corresponds to spins-up, dark 
to spin-down). B = 0.7J, d = 0.7J/spin and T = 0.45J. 

For a system with a sufficient large anisotropy it is 
energetically favorable to divide into parts with opposite 
directions of magnetization parallel to the easy axis in 
order to minimize the anisotropy energy barrier. This 
kind of reversal mechanism is called nucleation JI3] (see 
pj for a recent review). The simplest case of a reversal 
process driven by nucleation for a system of finite size 
is the growth of one single droplet starting somewhere 
at the boundary. Due to the growth of the droplet a 
domain wall will cross the system and the energy barrier 
which has to be overcome is caused by the domain wall 
energy. This is shown in Fig. || Here the domain wall is 
in the center of the system dividing the particle into two 
oppositely magnetized parts of equal size. 

The energy barrier AE n which has to be overcome dur- 
ing the reversal by nucleation of one single droplet in a 
three dimensional system with open boundary condition 
has been derived in || . It can be expanded with respect 
to MBR/a resulting in 



AE n w 2nR 2 a 
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(3) 



where a is the domain wall energy density. The corre- 
sponding lifetime of the metastable state for the single- 
droplet nucleation is then 



exp 



(AE r , 



V T 



(4) 



for temperatures T <C AE n . Eq. || has an interesting 
interpretation: for small fields B the energy barrier of 
a nucleation process is the energy of a flat domain wall 
in the center of the particle plus corrections which start 
linearly in B. In contrast to the Neel-Brown theory, here 
for vanishing magnetic field the energy barrier is propor- 
tional to the cross-sectional area of the particle rather 
than to its volume. 

Hence, comparing the two energy barriers for coherent 
rotation (Eq. g) and nucleation (Eq. |J) one can evalu- 
ate the critical particle size R c where the crossover from 
coherent rotation to nucleation sets in. For vanishing 
magnetic field this critical particle size is given by 



Rc ~ 2d' 



(5) 



For particles larger than R c reversal by nucleation has 
the lower energy barrier while the opposite is true for 
particles smaller than R c . 

Another possibility of a reversal process driven by nu- 
cleation is the growth of several droplets at the same time 
(multi-droplet nucleation). This reversal process occurs 
when the probability for the growth of large droplets is 
high enough and when the critical droplet size is smaller 
than the system size so that many droplets may occur in 
the system at the same time. This is the case for higher 
fields or larger temperatures jl5],[l6) . The lifetime for the 
multi-droplet nucleation in _D-dimensions is |L4fl 



exp 



/ AE„ 



\{D + 1)T 



(6) 



Hence, the exponent of the exponential function 
changes by a factor of 1/(D + 1) which is 1/4 in our case. 
In order to quantify the crossover from single droplet to 
multi-droplet nucleation we consider the case of infinite 
anisotropy in which case we can directly simulate an Ising 
system. We performed a standard MC simulation of the 
Ising Hamiltonian 



H= -J^SiSj - B^Si 



(7) 



with Si = ±1. 

Fig. H shows a corresponding spin configuration during 
the reversal. The Ising spins are represented as light 
(spin up) and dark (spin down) grey boxes. The reversal 
process is initiated by many droplets at the boundary of 
the system. Later all these droplets join each other so 
that the outer shell of the particle is reversed while the 
inner parts have still the initial magnetization direction. 
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FIG. 3. Snap-shot of a simulated spin system at the lifetime 
t for multi-droplet nucleation. Shown is an Ising system of 
size R — 8 spins. In order to show the interior of the system 
we have cut the sphere and show only one half of the system. 
B = 0.5J and T = 2.8J. 



droplet nucleation theory. Especially, the general result 
that the effective energy barriers for single- and multi- 
droplet nucleation are related by the D— dependent fac- 
tor is confirmed. A corresponding behavior concerning 
the field dependence of the the metastable lifetime was 
also observed recently 15 lq]. 

Acknowledgments: We thank M. Acharyya and 
K. D. Usadel for helpful discussions. The work 
was supported by the Deutsche Forschungsgemeinschaft 
through Sondcrforschungsbereich 166 and through the 
Graduiertenkolleg "Struktur und Dynamik heterogener 
Systeme" . 



le+06 
le+05 
le+04 
1000 
100 
10 

1 



1 


1 1 1 / 1 


1 1 w ■ 












B = 0.8 □ " 


1 


1 1 1 1 


B = 0.5 * 
i i i 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 

l/T [1/J] 



FIG. 4. Metastable lifetime r vs. l/T for two different 
magnetic fields. System size R — 4 spins. 

Fig. U shows the temperature dependence of the 
metastable lifetime following from simulations for two 
different fields. Each data point is an average over 100 
independent runs. The data show a linear behavior with 
different slopes for low and high temperatures, respec- 
tively. For low temperatures single-droplet nucleation 
occurs and we can compare the corresponding data with 
the theoretical prediction, Eqs. (|[ ||). The slopes of 
the straight lines shown are the energy barriers AE n ob- 
tained from Eq. || which have the theoretical values 40. 4 J 
(B = 0.5 J) and 23.5 J {B — 0.8 J) for the system size used 
here [R = 4 spins) and for a — 1.2J/bond || showing 
very good aggreement between the numerical data and 
the theoretical predictions. For higher temperatures we 
find a crossover to multi-droplet nucleation in which case 
theory predicts a reduction of the energy barrier by a 
factor of (D + 1) = 4, Eq. |[ The slopes of the straight 
lines shown are the energy barriers AE n / (D+l) obtained 
from theory showing again very good agreement with the 
simulated data. Note that the position of all lines are fit- 
ted since we do not know the prefactors in Eqs. [| and 
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To conclude we found very good agreement between 
the energy barriers obtained from the numerical data 
and the predictions obtained from both single- and multi- 
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